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Abstract. We calculate the bi-photon state generated by spontaneous parametric 
down conversion in a thin crystal and under collinear phase matching conditions using 
a pump consisting of any superposition of Laguerre-Gauss modes. The result has 
no restrictions on the angular or radial momenta or, in particular, on the width of 
the pump, signal and idler modes. We demonstrate the strong effect of the pump to 
signal/idler width ratio on the composition of the down-converted entangled fields. 
Knowledge of the pump to signal/idler width ratio is shown to be essential when 
calculating the maximally entangled states that can be produced using pumps with a 
complex spatial profile. 
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1. Introduction 

Quantum entanglement is one of the defining properties of quantum mechanics. It 
forms the basis for quantum information [K |2] and quantum computing [3j and is 
essential for applications such as entangled cryptographic systems [4J and quantum 
imaging [SJ [6l [71 |BJ • Much of the research in this area has focussed on qubit quantum 
entanglement, although entanglement has also been demonstrated between spatial 
modes carrying orbital angular momentum (0AM) ^ Uni HI]. Since these modes 
are defined within an infinite-dimensional, discrete Hilbert space there is considerable 
interest in their potential to generate multi-dimensional entangled states not only for 
enhancing the efficiency of current quantum protocols, but also for the realization of 
new, higher dimensional quantum channels. 

Spontaneous parametric down conversion (SPDC), the generation of two lower- 
frequency photons when a pump field interacts with a nonlinear crystal, has been shown 
to be a reliable source of photons which are entangled [12] and, in particular, in their 
0AM [9]. The spatial structure of the down-converted biphotons can be expressed as 
a superposition of Laguerre-Gauss modes of different amplitudes, with the width of the 
modal expansion relating to the amount of entanglement of the final state. In this 
paper, we calculate the exact analytical form of the biphotons for any Laguerre-Gauss 
pump, with no restrictions on the angular or radial momenta, i and respectively, or, 
in particular, on the width, of the pump, signal and idler modes. Our calculation 
demonstrates that 0AM is conserved. 

We use our result to calculate the coincidence amplitudes of the LG modes in the 
generated two-photon entangled state for a variety of pump beams. We demonstrate 
the the importance of the various beam sizes (the pump width and the choice of LG 
base for the signal and idler) in determining the state of the down-converted photon 
and the resulting £ distribution, also known as the quantum spiral bandwidth (SB) [13]. 
We find excellent agreement with previous analyses which considered more restricted 
conditions [131 E] • 

Knowledge of the exact form of the down-converted photons opens up the possibility 
of engineering the entangled modes, through an appropriate choice of pump beam, 
in order to produce biphotons which are, for example, maximally entangled in some 
subspace [I5]. We show how our result can be easily extended to pump beams which 
are complex superpositions of Laguerre-Gauss modes, including pump modes containing 
phase singularities. We again see excellent agreement with previous analysis [15j. 
Moreover, by allowing the beam widths to be free parameters, we show that the pump 
to signal/idler size plays a critical role in the engineering of the down-converted state, 
with states only being maximally entangled at particular width ratios. 
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2. Spontaneous parametric down conversion and spiral bandwidth 

We consider a thin nonlinear crystal (typically l-3mm), tuned for type-I, coUinear 
down- conversion, illuminated by a continuous-wave Laguerre-Gaussian pump beam 
propagating in the z-direction and confined within the crystal and assume perfect phase 
matching. This produces two highly-correlated, lower- frequency photons, commonly 
termed signal and idler. Since energy is conserved, ujp = ujg + uJi^ where the subscripts 
5, i refer to the pump, signal and idler, respectively. 

The down-converted biphoton state in this case is given by [I6l [T5] 

l^sPDc) = J rfr^$(r^)at(r^)aj(r^)|0,0) (1) 

where ^{r±) is the spatial distribution of the pump beam at the input face of the crystal, 
|0,0) is the vacuum state and aj(r^), aj(r^) are creation operators for the signal and 
idler modes, respectively. The integral is over the plane perpendicular to the axis of the 
pump beam. 

Since the photon pairs generated by SPDC are entangled in arbitrary superpositions 
of an infinite number of modes with 0AM they are best described by mode functions 
that are Laguerre-Gauss modes, LG^. At the beam waist {z = 0) the normalised LG 
modes are given, in cylindrical co-ordinates, by fl7j: 

Here w is the beam waist, i corresponds to the angular momentum, ih^ carried by the 
beam and describes the helical structure of the wave front around a wave front singu- 
larity and p + 1 describes the number of radial intensity maxima. 



(3) 



where P^^^^^- = 

Ps,Pi 



We can then write the biphoton state as: 

\i^spDc) = Y.Y.ct^:KPs■,e^,P: 

iPs iPi 

2 

is the probability of finding one photon in the signal mode \^s)Ps) 
and the other in the idler mode \^i^Pi) given a pump mode \^p^Pp). 

The coincidence amplitudes C^JJ^^ are calculated from the overlap integral 

'-'^ - {1P^M^SPDC) (4) 

27r (^OO r -I >k r i >k 

(5) 



Ps.Pi 



rZTT POO r- 1 * r 1 



Note that after substituting (Ej) into 

r27T 

) exp [i {£p — is — 



the integral over the azimuthal coordinate is 



/ 

Jo 



(6) 



from which we obtain the well-known conservation law for the orbital angular 
momentum: ip = is + U ^ IIH]. Substituting this into (Ej), and making use of the 
associated Laguerre polynomial in the form fT9] 

(p+KI)! 



i=0 



{p-l)\{\i\+l)\l\ 



X 



(7) 
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we obtain 



^nwl (1 + 72 + ^2^-^+1 



Pp Ps Pi (_0\k+i+j 2i '2j 

Esr^ sr^ \ Is H 
^—^ ^—^ /-I I 9 I 2\k+i+j 



[ae + k + i+j]l 



{Pp - mK\ + ky.kKPs - + tVAPi - jvm+m ' 

where = {\£p\ + l^^l + \£i\) /2 and 7^ = Wp/wg^ 7i = Wp/wi are the ratios of the pump 
width to the signal/idler widths. 

Equation ([8]) aUows us to calculate the exact analytical form of the down-converted 
photons produced using any Laguerre-Gaussian pump (or superposition thereof) with no 
restrictions on any of the beam parameters. Calculation of the coincidence amplitudes 
can also be used to calculate the number of 0AM modes participating in the down- 
converted state, otherwise known as the spiral bandwidth (SB). 



3. Coincidence amplitudes and spiral bandwidth for special cases 

For a given pump, equation ([8]) shows that the coincidence amplitudes have a strong 
dependence on the ratio of the pump width to the chosen widths of the LG bases for the 
signal and idler, 7^ and 7^. Experimentally, the maximum width of the pump is limited 
by the size of the down- conversion crystal, while the widths of the signal and idler beams 
are determined by the size of the single mode fibre used in the detection process. We 
can see this dependence more clearly if we consider some simplified examples. 



3.1. Gaussian pump, ip = Pp = 0. 

We first consider the simplest possible case! a Gaussian pump with £p = Pp = and 
biphotons with no radial dependence i.e. Ps = Pi = 0. Moreover, we assume that the 
signal and idler modes have the same width s.t. 7s = 7i = 7. From from conservation 
of 0AM, equation ([6]), we know that ig = —^i = ^ in this case and so we can reduce 
equation ([8]) to 



^0,0 — 



\ TTW'^ V 1 + 27^ / 

It is clear that the SB depends only on the ratio of the pump to signal/idler widths, 7. 
As in previous work, for example [131 K], this further reduces to 

CoVocflV (10) 



.3 

when all the beam widths are the same, i.e. 7 = 1. 
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The effect of 7 on the spiral bandwidth is shown in figure ([T]) where we have plotted 
the SB as a function of 7. It can be seen that as the signal/idler widths are reduced 
w.r.t. the pump width (i.e. as 7 increases) the SB increases but that the amplitudes of 
the participating modes decrease. In figure ([2]) we plot (in blue) the SBs for 7 = 1.0, 2.0 
and 3.0. These results are in excellent agreement with previous work, see, for example, 
figure (1) of [13] (shown in red in figure ([2])). Note that throughout this work although 
we only plot the SB from — 15 < ^ < 15, which are experimentally realisable values, we 
calculate the coincidence amplitudes for a much larger range of £ and use these results to 
"normalize" the coincidence probabilities. The effect of the different beam parameters 

3.0-| 

i.oJ 

0.5-lr ^ ^ ^ r 

-15 -8 8 15 

4 

Figure 1. Spiral bandwidth for a Gaussian pump, £p = pp = and fixed width 
k;^ = 1.0 as a function 7. 
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Figure 2. Spiral bandwidth for a Gaussian pump, ip = pp = O^Wp = 1.0. The signal 
and idler widths are equal and reduced s.t. 7 = 1.0, 2.0, 3.0 from left to right. Vertical 
(blue) lines are calculated using (j9j) and dashed (red) lines are calculated using eqn.s 

(10) & (11) of dg. 

with a Gaussian pump has recently been explored [20]. Note that the results are also 
valid for non-collinear geometries in the approximation of a thin crystal. 

3.2. Laguerre- Gaussian pump, £p > 0. 

We now consider what happens when the pump is a higher-order Laguerre-Gaussian 
mode, i.e. ip > 0. In this case conservation of 0AM requires £i = ip — £s- As before, we 
start with the simplest case: Pp^s,i = and assume that the signal and idler modes have 
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the same width. The probabihty amphtudes are then described by 



^0,0 



\ 



TTWt, 



Kp|!|4|!K«|!(l + 7H7f) 



where = + + |^^|) /2. This agrees with equation (8) of [14j when using 
normahzed LG modes and assuming that aU of the modes have the same width. 

To see how the 0AM of the pump affects the spiral bandwidth we first calculate 
the coincidence probabilities, Po,o ^~ % as a function of ig and 7 for different values of 
£p. As figure ([3]) shows, the SB is symmetric around the 0AM value of the pump. For 
small values of 7 the SB has a maximum at ig = ^p- As 7 is increased, however, the 



SB splits into two "wings" for > 0, with a minimum at 



and the amplitude 



increasing with ig- The physical explanation of this is that pump modes with ip > 
form a rings whose diameter increases with both the value of £p and the width of the 
beam. Maximum overlap between the pump and the signal/idler modes, and hence 
maximum coincidence amplitudes, therefore occurs for larger values of ig^^i- 

We next consider what happens if the widths of the signal and idler are different. 
We find that by varying the ratio of the signal to idler widths (in this case from 7^ = O.575 
to 2.O75) it is possible to fully suppress the SB for either negative or positive values of 
4, as shown in figure dH). 




Figure 3. Spiral bandwidth for a Laguerre-Gaussian pumps of fixed width Wp = 1.0 
and ip = 1, 2, 3, 4 as the signal/idler width is reduced s.t. 7^ = 7^ goes from 0.5 to 3.0. 



Until now we have considered only modes with no radial modes, i.e. Pp = Ps = 
Pi = 0. Higher-order LG beams are commonly produced using a spatial light modulator 
(SLM). Since these produce 0AM modes with a range of radial indices, [21j, we now 
look at the effect p has on the form of the spatial mode function of the entangled photons. 
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Figure 4. Spiral bandwidth for a Laguerre-Gaussian pump with ip = 2^Wp = 1.0 as 
the signal/idler width is reduced s.t. 7^ goes from 0.5 to 3.0. The signal and idler 
modes have different width ratios 7^/75 = 0.5, 0.75, 1.0, 1.5, 2.0. 

We find that the number of "wings" increases with the radial index of the pump as shown 
in figure ([5]). As before, altering the ratio of the signal to idler width suppresses one set 
of wings. Again these correspond to areas of maximum overlap between the pump and 
the signal/idler modes. 

If we allow first our signal/idler beams to have non-zero radial modes and then all 
beams to have non-zero radial modes the spiral bandwidth as a function of 7 becomes 
ever more complex. The effect of the radial number of the pump is shown in figure (|6j). 



3.3. Engineered entangled states 

For many quantum communication applications entangling systems in higher- 
dimensional states is important. Not only does the higher dimensionality imply a 
greater potential for applications in quantum information processing, but it has also 
been suggested that increasing the dimensionality of the entangled states of a system 
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Figure 5. Spiral bandwidth for a Laguerre-Gaussian pumps with £p = 2^Wp = 1.0 as 
the signal/idler width is reduced s.t. 7^ = 7^ goes from 0.5 to 3.0. The radial number 
of the pump is increased from = to 3. 
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Figure 6. Spiral bandwidth for a Laguerre-Gaussian pumps with ip = pp = 2^Wp = 
1.0 as the signal/idler width is reduced s.t. 7^ = 7^ goes from 0.5 to 3.0. The radial 
number of the signal and idler modes is increased from = = to 3. 
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can make its non-classical correlations more robust to the presence of noise and other 
detrimental effects [22[|^. Torres et al. demonstrated (theoretically) that arbitrary 
engineered entangled states in any rf-dimensional Hilbert space can be prepared using 
SPDC to translate the topological information contained in a pump beam into the 
amplitudes of the generated entangled quantum states p^. Controlling 0AM state 
superpositions in this way engenders the ability to produce and manipulate quantum 
states with an arbitrarily large number of dimensions. Indeed, Torres et al. illustrated 
this by calculating a number of maximally entangled states of different dimensions, 
produced using pump beams containing phase singularities (PS). Here we describe how 
our results can be easily extended to such pumps, which can be described by complex 
superpositions of LG modes [24j. Moreover, we demonstrate that the ratio of the pump 
to signal/idler widths must be taken into consideration in the calculation as it plays 
a critical role in the composition of the down-converted state: a state is generally 
maximally entangled for only one pump to signal/idler width ratio. 

We first calculate the coincidence amplitudes resulting from a pump that is a 
superposition of Laguerre-Gauss modes of equal width, and complex amplitudes 
a^. We can write this as 

^ a^LGll (p, 0) = a^LGl\ (p, 0) + a,LGl\ (p, 0) + . . . (11) 

n 

where n labels the modes in the superposition and |a^p = 1. The overlap integral 
([5]) then becomes 

r2'n roo ^ i * r i * 

^^:£=E / d4> pdpa^LGll{pA)[LGl{pA)\ [lG';XpA)\ 

J JO 

= ^^nC^pllp., (12) 
n 

where C^^^^'^^. are the coincidence amplitudes, calculated from ([8]), for each Laguerre- 
Gaussian component, \^p^Pp)•, of the pump field. 

In order to investigate a pump containing PSs we make use of its projection onto 
LG modes [25j. A field containing N PSs can be written as a superposition of LG 
modes, Y^LqCLiLGq^ with complex amplitudes 

a, = V^(-lf (^^y 'Vl^.b^.e (13) 
where bn are given by equation (11) of [15j or, equivalently, by 

''.. = ;^j^n(i+p.e'Mu (14) 

with m = [O^N] and Pi,(/>i being the radial and azimuthal positions of the PSs, 
respectively. 

Following [T5], we consider a pump containing six PSs at positions pi = 
0.65wp^ P2 1.85wp^ ps l.OGwp^ p4 O.bAwp^ p^ 1.531^;^, Pq 1.2Awp and 
0^ = i7r/3 for i = 1, 6. Substituting these values into f fT3l) and then into f fT2l) allows us 
to calculate the form of the pump, as shown on the left in figure ([7]), and the resultant 
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coincidence probabilities, vertical (blue) lines shown on the right in figure ([7j). We 
find excellent agreement with the results from eqn. (14) of [15j (dashed red line). We 
calculated the coincidence probabilities Pq'q ^ Pq^q ^ Pq^q , Pq'q and confirmed that these 
were equal and so this pump produces the maximally entangled qu-quart in the subspace 
5^4 = {|0, 0) , |1, 1) , |2, 2) , |3, 3)}. Note, however, that in this example we assumed that 




15 10 5 5 10 15 



4 

Figure 7. Pump containing six phase singularities, as described in [15 , and resulting 
probability amplitudes. The blue lines show the results found using equation (14) of 
[15] while the red line is the result found using equations (jS]) and (p!2]) . 

the pump, signal and idler all had the same width i.e. 7^ = 7^ = 1. Since we have seen 
that the composition of the down-converted state is strongly dependent on the ratio of 
the pump to signal and idler widths we extended these results by calculating the spiral 
bandwidth as a function of 7, keeping 7^ = 7^ for simplicity. As expected, we see a large 
variation in SB for different 7, as shown on the left in figure (|8]). This can been seen 
more clearly on the right of figure (|8]), where we plot the SB at 7 = 1.0, 2.0, 3.0. It is 
clear that by choosing the mode widths appropriately it is possible to tailor the spiral 
bandwidth in very different ways. 

Changing the relative beam sizes has a large effect on the composition of the 
entangled state with important implications in the "engineering" of the down-converted 
state: states which are maximally entangled for one value of 7 are generally not 
maximally entangled for any other values of 7 values. To illustrate this we consider 
the maximally entangled qu-quart produced by the six-PS pump above and calculate 
the coincidence probabilities Pq'q ^ Pq^q ^ Pq^q ^ Pq'q as a function of 7. From figure (|9j) it 
is clear that this state is only maximally entangled if 7 = 0.5 or 7 = 1.0. In order to 
obtain this state then, one has to ensure that these conditions are met. 

4. Conclusion 

We have calculated the exact analytical form of the biphotons produced from SPDC 
in a thin crystal and under collinear phase matching conditions. Our result has no 
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Figure 8. Spiral bandwidth for superposition pump for 7^ = 7^ = 0.5 to 3.0 (left) 
with slices at 7^ = 7^ = 1.0 (green, small dash), 2.0 (blue, large dash) and 3.0 (red, 
solid) shown on right. 
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Figure 9. Coincidence probabilities Pq'q (red, solid), Pq'q (blue, small dash), Pq'q 
(magenta, medium dash), and Pq'q (green, large dash) as a function of 7. These are 
equal, and so the state is maximally entangled, for only 2 values of 7: 7 = 0.5, 1.0. 

restrictions on the angular or radial momenta or, in particular, on the width of the 
pump, signal and idler modes and can be used with any pump that is a Laguerre-Gauss 
mode, or superposition thereof. We have shown that 0AM is conserved and found 
excellent agreement with previous analyses which considered more restricted conditions. 

We have used our result to calculate the spiral bandwidth for a variety of pumps and 
demonstrated the importance of the various beam parameters, in particular the pump 
width and the choice of LG base for the signal and idler, on the resulting coincidence 
amplitudes. 

In particular, we have demonstrated the importance of the beam widths when 
calculating engineered entangled states and shown that maximal entanglement is 
dependent both on the form of the pump and on the ratio of the beam widths. We 
suggest that this is an important consideration and may also offer another degree of 
freedom in quantum communications. 
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